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ABSTRACT 

We have performed infrared imaging of the jet of the quasar 3C 273 at wavelengths 3.6 and 5.8 /im with the 
Infrared Array Camera (IRAC) on the Spitzer Space Telescope. When combined with the radio, optical and 
X-ray measurements, the IRAC photometry of the X-ray-bright jet knots clearly shows that the optical emission 
is dominated by the high-energy emission component of the jet, not by the radio synchrotron component, as 
had been assumed to date. The high-energy component, represented by a power-law from the optical through 
X-ray, may be due to a second synchrotron component or to inverse Compton scattering of ambient photons. 
In the former case, we argue that the acceleration of protons exceeding energies of E p ~ 10 16 eV or possibly 
even to E p ~ 10 19 eV would be taking place in the jet knots of 3C 273 assuming that the acceleration time is 
proportional to the particle gyroradius. In contrast, the inverse Compton model, into which highly relativistic 
Doppler beaming has to be incorporated, requires very low-energy electrons of E e ~ 1 MeV in the jet knots. 
The present polarization data in the radio and optical would favor the former interpretation in the case of the 
3C 273 jet. Sensitive and detailed measurements of optical polarization are important to establish the radiation 
mechanism responsible for the high-energy emission. The present study offers new clues as to the controversial 
origin of the X-ray emission seen in many quasar jets. 

Subject headings: galaxies: jets — infrared: galaxies — quasars: individual(3C 273) — radiation mechanisms: 
non-thermal 



1. INTRODUCTION 

Large-scale jets extending to hundreds of kiloparsec dis- 
tances from the quasar nucleus often radiate their power 
predominantly in X-rays. This has been revea led by sur- 
veys with the Chandra X -ray Observatory ( Sambruna et al. 
120041 iMarshall et al.1120051) following its unexpected discov- 
ery of an X-ray jet in the quasar PKS 06 37-752 (z = 0.65: 
Chart as et all 120001 ISchwartz etaTJ 120001) . Despite exten- 
sive observational and theoretical work, the dominant X-ray 
emission mechanism operatin g in quasar jets remains unset- 
tled (e.g., lAtovan & Permed 120041 lHardcastlel 12006). Nei- 
ther naive single-component synchrotron nor synchrotron- 
self-Compton models can fully explain the spectral energy 
distributions (SEDs) at radio, optical and X-ray wavelengths. 
A widely discussed hypothesis for the strong X-ray emis- 
sion is relativistically enhanced inverse Compton (IC) scat- 
tering off the cosmic microwav e background (CMB) photons 
jTavecchio et all2000tlCelotti et all2001l) . In this model, the 
bulk velocity of the jet is assumed to be highly relativistic all 
the way to nearly megaparsec distances, with a Doppler factor 
of 8 ~ 10. Applying the beamed IC model to a large s ample 
of C/zaniira-detected quasar jets, Kataoka & Stawarzl (120051) 
found a range of 4 < 6 < 16 for ~ 90% of the sources. 

Alternative scenarios for producing X-ray emission at such 
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large distances from the central engine include a variety of 
"non-conventional" synchrotron models. The X-ray emis- 
sion could be explained by a spectral hump in a synchrotron 
spectrum as a result of reduced IC cooling in the Klein- 
Nishina regime provided that a Doppler factor of the jet 
is large (Derm er & Atovar] 120021) . Alternatively, a second 
synchrotron component responsible for the X-rays could be 
formed by turbulent accelera tion in the shear boundary layer 
( Stawarz & Ostrowski 2002) or by hypothetical high-energy 
neutral bea ms (neutrons and gamma-rays) from the c entral en- 
gine (Atovan & Dermer 2001; Neronov et al. 2002). Finally, 
synchrotron radiation may be produced by very-high-energy 
protons, E p > 10 18 eV ( Aharonian 2002). Discussion of these 
models is hampered mostly by the limited observational win- 
dows available so far. Apart from radio observations, only 
Hubble Space Telescope (HST) and Chandra have the capa- 
bility to resolve jets of powerful quasars. 

To explore the mid-infrared properties of the extended jet 
emission from powerful quasars and thereby to shed new 
light on the riddle of the C/ianiira-detectedjets, we performed 
Spitzer IRAC imaging observations of four powerful quasars 
during the first observing cycle. The first result of this pro- 
gram, based on the observation of quasar PKS 0637-752 
(Uchivamaet al. 2005), reported the first mid-infrared detec- 
tion of jet knots in a powerful quasar, and demonstrated that 
mid-infrared fluxes are of great importance in studying the 
SEDs, given the large gap bet ween the radio and optical wave- 
lengths. Following Georganopoulos et al. (2005), we also 
showed that in the context of the beamed IC model, observa- 
tions with the Spitzer IRAC c an constrain the matter content 
of jets dUchivama et all2005l) . 

In this paper, we present Spitzer IRAC results for the well- 
known jet in the nearby bright quasar 3C 273. We combine 
Spitzer infrared measurements with available data from the 
Very Large Array {VIA), HST, and Chandra, forming an un- 
precedented data set to study the physics in extragalactic jets, 
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FIG. 1 . — Multifrequency images of the 3C 273 jet normalized to the the maximum pixel, (a) Spitzer IRAC (PSF-subtracted) image at 3.6 fim; (b) HST image 
at 620 nm with VLA (at 2.0 cm) contours overlaid [adapted from Jester et al. 1 2005)]; (c) Chandra image in the 0.4-6 keV interval. The origin of the coordinate 
system is at the 3C 273 core and the horizontal axis, parallel to the ridge line of the jet, is set at position angle 222° (the quasar is to the left at x = 0). 



not easily attainable for any of the other C/zawc/ra-detected 
jets found in more distant quasars 9 . With these data we show 
that the multiwavelength SED can be decomposed into two 
distinct nonthermal components. Based on our finding, we 
consider plausible scenarios for the origin of the X-ray emis- 
sion, and derive quantitative constraints on the physical state 
of the jet knots. 

This paper is organized as follows. In |2]we summarize 
previous results on the 3C 273 jet. In 5j3] we present analysis 
of the new Spitzer IRAC data and describe the radio, opti- 
cal, and X-ray data for the same emission regions. The mul- 
tiwavelength spectral energy distributions are presented and 
discussed in ^ Concluding remarks a re given i n §5\ 

The redshift of 3C 273 is z = 0.158 (Schmidt 1963), so we 
adopt a luminosity distance of Di = 749 Mpc (at this distance 
1" corresponds to a linear scale of 2.7 kpc) for a ACDM cos- 
molo gy with tt m = 0.27, p A = 0.73, and H = 71 km s" 1 Mpc" 1 
(e.g.. lSpergel et alJ20()3). We define the spectral index, a, in 
a conventional way according to the relation f u oc v~ a , where 
f v is the flux density at frequency v. 

2. A BRIEF SUMMARY OF THE 3C 273 JET 

The presence of the kiloparsec scale jet of 3C 273 was 
recognized already in the famous quasar-discovery paper of 

9 A list of C/MMt/ra-detected extragalactic jets is found at 
http://hea-www.harvard.edu/XJET/ 



Schmid] dl963h. The iet of 3C 273 has bee n extensively 
studied (see Courvoisier 1998; Stawarz||2004t for r eviews): 
in the radio JDavis et all 1 19851 [Con way et alJll993l). in the 
near-infrared {M cLeod & Rieka l 1 994T] Neumann et alJI 19971 
|Jester etay |2005 j). i n the optical jRo ser & Meisenheimer 
1 199 It iBahcall et alJ H995t Hester et al] 1200 lL 120051) and 
in the X-ravs llHarris & Sternl 119871: iRoser et alJ l20"o"(l 
iMarshall et"afll200 It ISambruna et al JI200 ll) . Fi gure fTl shows 
the jet at various wavelengths (see 33.2l for details). Radio and 
optical emission from the large-scale jet are almost certainly 
synchrotron in nature given their measured polarizations and 
the coincidence in orientation. However, the emission mech- 
anism giving rise to the X-rays has been controversial (see 
below). 

The one-sided radio jet appears continuous all the way from 
the nucleus to t he "head" region inc reasing in brightness to- 
ward the head (Conwav et al. 1993). The head region is lo- 
cated ~ 23" from the nucleus or a projected distance of 62 
kpc. The jet shows a sideways oscillation or "wiggle" with 
respect to the ridge line of the jet. The nonthermal radio 
spectrum is characterized by spectral index a la dio — 0.7-0.8; 
the spectral index shows little variation along the jet. At the 
southeastern side the jet is accompanied by an extended but 
narr ow "lobe" that sh ows a quite steep spectrum of a ra dio ~ 
1.5 ( D avis et aiTT985). The faint lobe may be backflow from 
the hot spots or may represent outer boundary layers. 
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Opti cal imaging of the large-scale jet with HST 
( Bahc aTTet al.l 119951) revealed a structured morphology 
that is characterized by successive bright knots interspersed 
by regions of weak emission. The optical jet, extending from 
knot A to H2 10 , is coincident with the radio features, but the 
brightness profile in the optical differs from that in the radio 
frequencies . The optical spectra are hard er in regions close 
to the core (Roser & Me isenheimerl 1 99"fl) . in such a way that 
spectral index declines from aw ~ 0.6 at knot A to a opt ~ 1 .5 
at D2/H3 (Jeste r et al.ll2001l) . It has been demonstrated that 
the flux at 300 nm (near ultraviolet) s hows an excess ov er the 
extrapolation from lower frequencies (Jeste r et al.l2002l) . 

Both the radio and optical jet emission are linearly polar- 
ized. From x = 15" to 19" (where x denotes the angular 
distance from the core) the degree of radio polarization is 
10-20%, and the average is consistent with the optical value 
(Roser & Meisenhei merl 1 99 ltlRoser et alJl9 96). This agree- 
ment argues that both the optical and radio emission have a 
common (synchrotron) origin. The inferred magnetic field is 
longitudinal along the ridge line of the jet out to x ~ 21", be- 
yond which the sudden transition of the magnetic field config- 
uration from longitudinal to tra nsverse is observed in the po- 
larizati on data in both the rad i o dCo nwav et all 19931) and op - 
tical dScarrott & RolpU[l9"89l IRoser & M eisenhei me3ll991l) . 
Based on the distinctive polarization pattern and the radia- 
tion spectrum, the head region is widely considered to harbor 
terminal shocks and, as such, is physically different from its 
upstream knots. 

Excess X-ray emission associated with the jet has been 
reported based on observations made with EIN STEIN 
( lHarris & SterrJ 19871) and ROSAT (IRoser et all2000l) . Chan- 
dra observations have confirmed the X-ray emission along 
the jet and allowed detailed study of its spatial and sp ectral 
properties iSambruna et alJl200lt iMarshall et al.ll2001l) . X- 
ray flux is detected along the whole length of the optical jet 11 , 
from x ~ 12" (knot A) to ~ 21". Regions closer to the core, 
specifically knots A and B2, are particularly bright in the X- 
ray. The X-ray spectrum shows a power-law continuum, with 
spectr al index q x ~ 0.6-0 .9, indicative of its nonthermal na- 
ture llMarshall et al]l200"ll) . For knot A, a synchrotron-self- 
Compton (SSC) model (with an equipartition magnetic field 
and without Doppler beaming) predicts the SSC X-ray flux 
low er by three orders of magnitude than the observed flux 
/see IRoser et al.ll2000l ISamhSma et al JI200U IMarshall et al J 
120011) . 

It has been claimed that the innermost knot A has a straight 
power law continuum extending from the radio to X-ray, 
while the subsequent knots are described by a power law 
with a cut-off around 10 1 4 Hz plus a separate X-ray com- 
ponent (IRoser etal J 120001 IMarshall et alJ 120011) . If so, the 
radio through X-ray emission in knot A may be ascribed to 
synchrotron emission from the same electrons as the radio 
emission; the X-rays in the other kn ots are then produced in 
some other way. On the other hand, Sambr una et alJ (2001) 
argued that knot A as well as subsequent knots has a spec- 
tral break at optical wavelengths and suggested the SED is 
comprised of two components, a synchrotron part (radio to 
optical) and a beamed IC part (X-ray an d beyond). F inally, 
based on deep VLA and HST observations. lJester et alJ ( 120021) 




10 In this paper, following Jester et al. 1 2005), we adopt a nomenclature of 
knots as depicted in Fig.0x Note that there are many versions of the knot 
labeling in the literature. 

1 1 The presence of fain t jet emissio n closer to t he core (x < 10") is k nown 
in both the optical i Martel et al. 2003) and X-ray i Marshall e t all2001l) . 
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FIG. 2. — Spitzer IRAC image of 3C 273 and its jet at 3.6 fim in units of 
MJy sr . North is up and east is to the left. The origin of the coordinates is 
at the quasar core. The jet emission is clearly visible in the box drawn in the 
lower right quadrant. 



found a significant flattening in the near-infrared to ultravio- 
let spectra and argued that the radio-optical emission cannot 
be ascribed to single-population synchrotron radi ation requir- 
ing an additional high-energy component (see also lJester et alJ 
120051) . Below we construct the SEDs of the jet knots based on 
the extensive VLA, Spitzer, HST, and Chandra data, which al- 
low us to elucidate the controversy in interpreting the shapes 
of the broadband SEDs. We then discuss possible origins for 
the high-energy emission, X-rays in particular. 

3. OBSERVATIONS AND ANALYSIS 

3.1. Spitzer IRAC Observation 

The infrared observations of 3C 273 with Spitzer IRAC 
were carried out on 2005 June 10 as part of our GO-1 pro- 
gram (Spitzer Program ID 3586). IRAC is equipped with 
a four-channel camera, InSb arrays at wavelengths 3.6 and 
4.5 /im and Si:As arrays at 5.8 and 8.0 iim, eac h of which 
has a field of view of 5.'2 x 5.'2 ( iFazio et alf2 004). The InSb 
arrays operate at ~ 15 K while the SiAs arrays operate at 
~6K. Only the pair of 3.6 and 5.8 /im arrays, observing the 
same sky simultaneously, was chosen for our observation, to 
obtain longer exposures in one pair of bandpasses as opposed 
to two pairs with half the exposure time. The 3.6/5.8 pair 
was chosen because of its better spatial resolution and sensi- 
tivity. The pixel size in both arrays is ~ 1"22. The point- 
spread functions (PSFs) are characterized by widths of 1 "66 
and 1 "88 (FWHM) for the 3.6 and 5.8 fim bands, respectively 
( IFazio et al J 12004) . The flux measurement of p oint sources 
with IRAC is calibrated to an accuracy of 2% ( Rea ch et alJ 
2005). 

For each IRAC band, we obtained 26 well-dithered frames, 
each with a 30 s integration time. The pipeline process (ver- 
sion S 12.0.2) at the Spitzer Science Center (SSC) yielded 
26 calibrated images. These calibrated frames were com- 
bined into a mosaic image wi th a pixel size of 0."2 us - 
ing the SSC package MOPEX (Mako voz & Marleaull2"o"05l) . 
which removes spurious sources such as cosmic rays and 
moving objects based on inter-frame comparisons. The 
dither ed frames were co mbined using the "drizzling" tech- 
nique (Fruchter & Hook 2002). Finally, we made use of the 
astrometric positions of three field optical sources that have 
infrared counterparts to accurately (< 0."2) locate the position 
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of 3C 273 and its jet in the IRAC image. The 3.6 /mi mosaic 
image is shown in Fig. [21 Note that although the quasar core 
is extremely bright, the jet emission clearly stands out to the 
southwest. 

As is evident from Fig. |2] the infrared emission of the jet 
is contaminated by the PSF wings of the bright quasar core. 
Therefore, measuring the infrared fluxes of the jet compo- 
nents requires careful removal of the wings of the core un- 
derneath the jet emission. Unlike the case of PKS 0637—752 
dUchivama et alJ2005l) . the point-spread function provided by 
the SSC is not applicable for the much longer jet of 3C 273. 
Instead we made use of deep IRAC images of a very bright 
star in the Extended Chandra Deep Field South (E-CDFS) to 
generate a PSF template that is bright enough out to a dis- 
tance of ~ 20" away from the core. In fact, the template star 
is brighter than 3C 273 itself by a factor of ~ 13 in the 3.6 /im 
band and by a factor of ~ 4 in the 5.8 /mi band, respectively. 
Due to the saturation in both the 3C 273 and E-CDFS images, 
and possible contributions from the host galaxy of 3C 273, the 
PSF subtraction can only be carried out reliably for x > 10", 
but fortunately it covers the entire length of the main body 
of the optical and X-ray jet. The PSF subtraction procedure 
introduces an additional systematic error in the infrared flux 
determination. Based on the small level of intensity due to a 
residual in PSF subtraction in the vicinity of the jet, we have 
estimated the systematic errors to be 0.05 MJy sr -1 for the 
3.6 /im channel and 0.3 MJy sr -1 for the 5.8 /im channel, 
which are added in quadrature to the noise of each pixel when 
photometry is performed. The jet is so bright that photometry 
errors due to noise become unimportant with respect to these 
systematic errors. 

The core emission of 3C 273 itself is saturated in both the 
3.6 /im and 5.8 /im science frames. An additional frame was 
taken in a high-dynamic-range mode with a short integration 
time (1 s), such that 3C 273 itself is unsaturated; we use this 
frame to estimate the mid-infrared flux of the quasar core. 
We obtain the flux density of the 3C 273 core to be ~ 170 
mJy at 3.6 /im and ~ 230 mJy at 5.8 /im, respectively. At 
the relevant wavelengths, the spectrum of 3C 273 is known 
to have a steady small bump at ~ 3.5 um s uperposed on a 
variable power-law continuum (Robson et al. 1986). At the L 
band (3.45 um). lTurler et alJ (11999) list the mean flux density 
of f v = 160.5 mJy with a dispersion (due to time variabil- 
ity) of cr„ = 33.5 mJy, with which the IRAC flux is consis- 
tent. The spectroscopic observation of 3C 273 on 2004 Jan- 
uary 6 made with the Infrared Spectrograph (IRS) on Spitzer 
gives the flux density of /„ ~ 190 mJy at 5.8 /im llHao et alJ 
2005). The power-law continuum in June 2005 (our obser- 
vation) seems slightly brighter than January 2004. It should 
be noted, however, that the flux estimate from a single short 
frame suffers from larger systematic errors as compared to the 
science frames. 

3.2. Multiwavelength Jet Images 

Figure^ shows the IRAC image at 3.6 /im of the 3C 273 jet 
after subtraction of the PSF wings of the quasar core, where a 
line in position angle 222° from the core is adopted as a refer- 
ence horizontal axis. The 5.8 /im image is essentially similar 
to the 3.6 /tm image but with slightly worse resolution. In 
Fig. ^ we a ls° show high-resolution radio (at a wavelength 
of 2.0 cm obtained with VIA), optical (at 620 nm with HST), 
and X-ray (with Chandra) images. The VIA and HS T im- 
ages are same as those presented in lJester et alJ {2005), both 
of which have an effective resolution (FWHM) of 0."3. (The 
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FIG. 3. — Intensity profiles of the 3C 273 jet along the ridge line integrated 
over the vertical dimension from —1" to 1" (a) at 3.6 fim and (b) at 5.8 fim 
obtained with Spitzer IRAC, together with the best-fit models comprised of 9 
knot features. 



resolution in the HST image is artificially degraded from its 
native resolution to facilitate a meaningful comparison with 
other bands. The strongest radio spot H2 is truncated in the 
VIA contour map.) The Chandra image is based on new and 
deep observations with the ACIS-S detectors (PI: Jester). We 
des cribe the data analysis of the Chandra observations later 
in 33.51 [see also Jester et al. (2006, submitted) for comple- 
mentary analyses]. Though the X-ray image presented here is 
essentially identical to Fig. 1 of Marsha ll et alJ (12001). which 
was derived from the early calibration observations with the 
grating detectors, the spectral parameters are updated taking 
into account calibration improvement. 

The IRAC infrared emission traces almost perfectly the 
main body of the optical jet; the onset of the infrared jet coin- 
cides with the innermost optical knot A centered at x ~ 13", 
and the infrared jet starts fading around x ~ 21" in the head 
region. A wiggling shape with a transverse amplitude of 
Ay ~ 0"2 that is evident in the optical can barely be seen in 
the 3.6 /im image. The infrared jet brightness shares an over- 
all spatial pattern with the radio map in a sense that the bright- 
ness of the jet increases toward the outer part almost monoton- 
ically. However, the radio intensity peaks at H2 where the in- 
frared jet has already started getting fainter. The overall jet in 
the IRAC bands looks more similar to t hat at the near-infrared 
wavelength (see Neuma nn et alJ[l997t for a K' band image). 
In contrast, in the optical at 620 nm each knot is more or 
less of similar brightness, suggesting a spectral change from 
longer wavelengths. 

To further investigate the jet emission quantitatively, partic- 
ularly in terms of multiwavelength analysis, we derive the flux 
densities of individual knot components, at radio, infrared, 
optical, and X-ray wavelengths. The results of the photome- 
try are summarized in Tabled and presented in a conventional 
vf v form in Figs.|5]and|S] In the subsequent sections, we de- 
scribe the methods of our measurements in each band before 
proceeding with multiwavelength analysis in |4] 

3.3. Infrared Photometry 

The angular separation of adjacent knots (typically ~ 1") 
is comparable to the width of the IRAC PSF, which allows 
us to derive knot fluxes nearly individually even though the 
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TABLE 1 

Flux Densities of 3C 273 Jet Knots 



Flux Density, f \ 



Frequency 
i/(Hz) 


A 


Bl 


B2/B3 


CI 


C2 


Dl 


D2/H3 


H2 


HI 


VLA 1 (mjy): 




















8.33 X 10 9 


52 


42 


77 


70 


93 


250 


500 


750 


330 


1.50 X 10 10 


33 


27 


54 


48 


63 


160 


310 


440 


200 


2.31 X 10 10 


24 


20 


38 


36 


45 


110 


220 


280 


130 


Spitzer b (/xjy): 




















5.23x10° 45±10 


35±12 


49±11 


98±11 


89±12 


154 


161 


87±12 


28±11 


8.45 x 10 13 


27 


15±2.6 


39 


36 


46 


80 


140 


41 


10±2.6 


HSr OuJy): 




















1.87 x 10 14 


7.6 


4.4 


11 


9.0 


11 


21 


30 


5.4 


1.5 


4.85 x 10 14 


3.6 


1.7 


4.0 


2.6 


2.5 


3.6 


6.2 


1.1 




1.00 x 10 15 


2.4 


0.93 


2.0 


1.1 


0.90 


1.2 


2.1 


0.40 




Chandra 6 (nJy) 




















2.14 X 10 17 


46 


11 


26 


5.8 


5.1 


6.1 


7.6 






3.77 X 10 17 


28 


6.2 


14 


2.7 


2.6 


3.7 


3.9 






1.05 x 10 18 


12 


2.7±0.3 


5.1 


0.9±0.2 


1.1±0.2 


0.9±0.2 


1.4±0.3 







NOTE. — The 90% errors are shown only when exceeding 10%. Otherwise a formal error of 10% is assigned 
to conservatively account for the systematic errors (see the text). The optical fluxes are corrected for Galactic 
extinction using Ay = 0.07. The X-ray fluxes are corrected for Galactic absorption using Nu = 1 .7 x 10 20 cm -2 . 

"Nominal wavelengths are A = 3.6, 2.0 and 1.3 cm. 
b Nominal wavelengths are A = 5.73 and 3.55 fim. 

c Nominal wavelengths are A = 1.6, 0.62, and 0.3 fim. Knot HI is detected in the 1.6 fim band, but not at shorter 
wavelengths Jester et al. 2005). 
d Nominal energies are e = 0.89, 1.56, and 4.35 keV. 



infrared images do not show discernible knots. We have mod- 
eled the two-dimensional brightness distribution of the whole 
length of the infrared jet with a series of PSFs of variable in- 
tensity centered on the positions of the known optical knots. 
We constructed the PSF based on nearby template stars. To 
account for possible slight variations in the widths of the PSFs 
by ~ 10%, we adjusted each PSF width so that the PSF shape 
fits well the jet in the y-direction. As shown in Fig. \l]p, 1 1 
knots are predefined from the optical image 12 . Among them, 
B2 and B3 were combined into a single component since the 
two knots are too close (separated only by 0."5) to be fit indi- 
vidually. Similarly, D2 and H3 were united. Elongated PSFs 
were applied to the combined knots so as to account for the 
multiplicity. Whereas the optical spectra of D2 and H3 are 
similar to one an other, the optical slopes of B2 and B3 are 
slightly different tTester et al]l2005l) : one must therefore be 
careful when interpreting the summed spectra of B2/B3 . Also , 
we note that an infrared-bright star, St 2 in LTesteretalJ 12005). 
would contribute about 5% of the flux density of knot Dl. To 
summarize, we fit the jet images at 3.6 and 5.8 /xm with a 
series of 9 point-like components at fixed positions, after con- 
volving with the IRAC PSF. 

The PSF fit gives the flux density for each of the 9 com- 
ponents (Table 0. Note that in determining the photometric 
errors we have included the uncertainties imposed by back- 
ground (the wings of the core and the host galaxy) subtrac- 
tion. Also, photometry is subject to systematic errors due to 

12 We use a term "knot" for labeling but faint inter-knot emission as well 
as the so-called cocoon are included in our photometry. 



A Bl CI C2 Dl H2 m 

B2/B3 D2/m 



FIG. 4. — Apertures used for radio and optical photometry, superposed on 
the HST jet image (same as Fig.fTt). 



different spectral shapes, as the standard IRAC calibration is 
strictly valid for the specific case of /„ oc v~ l . Though we do 
not correct for this effect, the flux errors in this regard are less 
than 1 % for the spectral slopes of interest. 

Figure [3] displays the one-dimensional profile of the in- 
frared jet emission along the ridge line of the jet. The mea- 
sured brightness profile and the best-fit 9-component model 
are shown, integrated over the vertical dimension from -1" to 
1", together with its decomposition into individual PSFs. 

3.4. Radio and Optical Photometry 

We also performed photometry at three radio wavelengths 
{VIA: 3.6, 2.0, and 1.3 cm) and three optical wavelengths 
(HST: 1.6 fxm, 620 nm, and 300 nm). T he data sets used fo r 
the photometry are identical to those in iJester et alJ 1 2005). 
We re -analyzed the photometric images because Jest er et alJ 
d2005l) performed photometry with a 0"3-aperture while we 
need knot-by-knot photometry to be directly compared with 
the infrared analysis. Unlike the infrared photometry de- 
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scribed above, a PSF-fitting procedure is not mandatory since 
the knot separation of ~ 1" is large compared to the spatial 
resolution. Therefore, defining regions as depicted in Fig.[4l 
we extracted the corresponding radio and optical fluxes in- 
side each aperture region. The optical fluxes are corrected for 
Galactic extinction (Ay = 0.07); correction factors used here 
are 1.01 for 1.6 /im, 1.05 for 620 nm, and 1.12 for 300 nm, 
respectively. Photometric errors are small (< 5%) in all cases 
for both the radio and optical measurements. Possible system- 
atic errors such as those associated with deconvolution in the 
interferometric radio data or flux contamination from adjacent 
knots, though difficult to quantify, would be several percents. 

The morphologies of the knots in the jet are similar at all 
wavelengths; only the relative intensities of the knots change 
with frequency. The only exception to this is knot Bl, in 
which the radio emission does not have a corresponding fea- 
ture to the optical knot. This issue should be kept in mind 
when interpreting the spectrum of knot B 1 . 

3.5. X-ray Photometry 

Data used for the X-ray photometry presented here are 
based on four Chandra ACIS-S observations performed be- 
tween November 2003 and July 2004 with Observation IDs 
(Obs IDs) from 4876 to 4879. The aimpoint was placed at 
the S3 CCD. The basic data reduction was done by using the 
CIAO software version 3.2.2 and the CALDB 3.1.0 with the 
standard screening criteria applied. Periods of background 
flaring often caused by the increased flux of charged particles 
in orbit have been excluded from each data set (e.g., a ~ 6 ks 
flare was found in the data of Obs ID 4879). After the screen- 
ing processes a total exposure of 136 ks remains. We define 
three energy bands as 0.6-1 keV (soft), 1-2 keV (medium), 
and 3-6 keV (hard). The nominal energy e n orm for each band 
is then the mean energy weighted by the effective area: e n orm = 
0.89 keV, 1.56 keV, and 4.35 keV, for the soft, medium, and 
hard band, respectively. 

We have constructed X-ray flux images in the three bands 
with a pixel size of 0"2 by dividing a counts map by an 
exposure map. These flux images, which are actually the 
weighted mean of four observation sets, have units of photons 
cm -2 s" 1 keV" 1 per pixel and allow us to perform X-ray pho- 
tometry. Calculations of the exposure maps were performed 
with the CIAO tool mkexpmap assuming an incident spectrum 
of a power law with spectral index a = 1 , correcting for the 
contamination buildup in the ACIS detectors. The results are 
insensitive to a since the corrections for different indices are 
minimal at the effective-area-weighted mean energy. Then we 
have measured the flux densities at e norm by integrating over 
relevant apertures and multiplying each image by e n nrm- A 
similar but different procedure was adopted by Harr is et alJ 
(2004) to obtain the flux densities of the jet in 3C 120. Com- 
pared to the radio and optical cases (Fig. |4j, we elongated 
the object apertures in the transverse y-direction (\y\ < 1") 
to accumulate photons such that the portion of photons ly- 
ing outside the aperture should be < 10%. We estimated the 
background using regions outside the jet in the vicinity of the 
extraction apertures. The background amounts only to < 2% 
of the source flux for knot A even in the hard band where it 
has the largest contribution. Finally, to c orrect for the inter - 
stellar absorption of A^h = 1 .7 x 10 2() cm" 2 (lAlbert et alJ 19931) . 
we multiplied the measured fluxes by a factor of 1 .08 and 1 .02 
in the soft and medium energy bands, respectively. The flux 
densities of the knot regions determined in this way are given 
in Table [2 Most of the measurements have small statistical 



uncertainties less than 10%. The jet emission downstream of 
knot B2 does not show well discernible peaks but emission 
features associated with knots CI, Dl, and D2/H3 can be rec- 
ognized in the X-ray image, which justifies to some extent our 
aperture photometry for the downstream knots. 

In addition to the photometry-based spectra, we have per- 
formed the more usual spectral analysis as a consistency 
check. The X-ray spectra in the 0.5-6 keV range of the four 
observations were simultaneously fitted by a power law with a 
fixed absorption column A^h = 1.7 x 10 20 cm" 2 . The resultant 
spectral index for knot A is found to be a = 0.84 ± 0.04. Also, 
the index derived for a combined downstream region from 
knot CI to D2/H3 is a = 1.06 ±0.07. The photometry-based 
indices, on the other hand, are a ~ 0.8 and a ^ 1.0 for knot 
A and the combined downstream, respectively. Both meth- 
ods are thus consistent with one another. It should be noted 
that we obtained somewhat steeper spectra, by Act ~ 0.2 , 
compared to the results presented by Marshall et al. (2001). 
The difference arises most likely as a result of absorption by 
the molecular contaminants in the ACIS detectors. The pres- 
ence of such c ontaminants was not known at the time of the 
Marshall et al. ( 2001 ) analysis and consequently their spectral 
analyses inevitably led to a systematically lower a. Similarly, 
Perlman & Wilson (2005) reanalyzed early Chandra data on 
the M 87 jet and found that contamination of the ACIS de- 
tector can account for the flatter X-ray spect ra than ex pected 
based on the optical slope as reported by IWilson & Yang l 
(2002, see their erratum). 

4. RESULTS AND DISCUSSION 
4.1. SED of Jet Knots: Two-Component Nature 

In Fig. |3 we show the SEDs constructed for knots A, Bl, 
CI, C2, and Dl. Since the interpretation of multiple knots 
(B2/B3 and D2/H3) could be misleading, the "summed" spec- 
tra for B2/B3 and D2/H3 have been omitted in Fig. |5] and in 
the spectral modeling described below. The overall spectral 
appearance of the inner knots (A and Bl) has strikingly dif- 
ferent characteristics as compared to the outer knots (CI, C2, 
and D 1 ). First, the inner knots radiate strongly in the X-ray (or 
even beyond) with hard spectra while the outer knots radiate 
predominantly in the mid-infrared. Second, and most impor- 
tantly, the SEDs of the inner knots have a concave shape from 
the mid-infrared to optical band, indicating that two spectral 
components cross over at optical wavelengths. Figure[6]gives 
a close-up version of Fig. [5] highlighting the mid-infrared to 
near-ultraviolet part of the spectrum. We emphasize here that 
the IRAC fluxes fill a critical region in the SEDs, making 
it possible to identify unambiguously the spectral structures 
present in the infrared to optical wavelengths. 

In the SEDs of the inner knots, we identify two spectral 
components, namely (1) the low-energy synchrotron spec- 
trum extending from radio to infrared with a spectral cutoff 
at ~ 5 x 10 13 Hz, and (2) the high-energy component arising 
in the optical and smoothly connecting to the X-ray flux. Ten- 
tatively, we suggest that the high-energy component is likely 
of synchrotron origin, because the optical polarization is con- 
siste nt wi th t he ra dio polarization in degree and orientation. 
In § 34.31 and 14.41 we discuss the origin of the high-energy 
spectral component and its implications to particle acceler- 
ation in the jet. In explaining the broadband spectrum of knot 
A, a single power-law continuum all the way from radio to 
X-ray has been prop osed in prior studies ( Rtiser et al .1 120001 
Marsh all et alJ l2001). but this picture does not correctly de- 
scribe the present, more extensive data. 
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FIG. 5. — Broadband SEDs of the 3C 273 large-cale jet. (a) Inner knots A (filled circles) and Bl (open circles); (b) outer knots CI (diamonds), C2 (open 
circles), and Dl (filled circles). The photometric data are given in TablelTI The curves represent the two-component phenomenological model, namely equation 
ffTl : f v = K[V~ ai exp [- (f/v c ) a ] +K2^° 2 with parameters in Tablel2l Note that 1 Jy Hz = 10~ erg cm -2 s -1 . PigurelSlzooms in the mid-infrared to ultraviolet 
part, where the photometric points are relatively cramped. 




In the outer knots, the SEDs also show the two components 
but with different relative strengths. The low-energy part 
becomes more dominant with increasing distance from the 
quasar core. At first glance the entire infrared-optical spec- 
trum may be attributed to the falling part of synchrotron emis- 
sion. However, careful inspection shows a slight spectral- 
flattening in the optical and near-ultraviolet, which suggests 
a contribution from the second component responsible for th e 
X-ray emission, as already proposed by Uester et all fc002l) . 
This idea is further strengthened by the observation that the 



radio brightening in Dl (relative to C1/C2) accompanies in- 
frared brightening while the near-ultraviolet and X-ray fluxes 
remain at similar levels (see Fig.[6j. 

To visualize the spectral transition from the inner to the 
outer knots along the jet, we present a three-color image of 
the jet in Fig. \J\ based on the data taken with Spitzer, HST, 
and Chandra, with the VIA contours from Fig. ^ overlaid. 
The colors are coded as follows: Spitzer 3.6 /jm {red), HST 
/o.3Mm _ 0.05 x fi.6fj.rn (green), Chandra 0.4-6 keV (blue). The 
Spitzer image is a "deconvolved" IRAC image, which repre- 
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sen ts th e results of our infrared photometry via PSF-fitting 
in SI3.3I Specifically, the best-fit image obtained with two- 
dimensional fitting is artificially shrunk to match the effec- 
tive resolution of the HST image. The HST image is an "UV 
excess" m ap [determined as fm„,m -0.05 x fit^m based on 
the data in Jeste r et all (1200 5)1 representing the dominance of 
near-ultraviolet light over the near-infrared. It is clear that the 
inner knots are luminous in both the UV excess and X-rays 
(the high-energy component), while the outer knots are bright 
in the mid-infrared (the low-energy component). 



TABLE 2 

Phenomenological Model Parameters of 3C 273 Jet Knots 



Parameter 


A 


Bl 


CI 


C2 


Dl 


H2 


HI 


at 


0.76 


0.73 


0.63 


0.68 


0.75 


0.9 


0.9 


v c (Hz) 


5 x 10 13 


5 X 10 13 


4x 10 13 


5 X 10 13 


6 X 10 l: 


5 x 10 13 


4x 10 13 


a 


0.8 


0.8 


0.8 


0.8 


0.8 


0.8 


0.8 


OL2 


0.7 


0.8 


1.0 


1.0 


1.0 







NOTE. — Equation Q, /„ = K\v a < exp [- (i//i/ c ) a ] + kjv 
adopted to characterize the SEDs. 



We describe the radio-to-X-ray spectra with the following 
double power-law models (one with an exponential cutoff): 



U =f\v + fiv = K\v Q > exp [- {yl v c ) a ] + n 2 v 



(1) 



The first term of the right-hand-side accounts for the low- 
energy spectrum and the second pure power-law describes the 
high-energy part. (For the second component, f 2v = n 2 v~ ai , 
we set an artificial low-energy turn-over below 10 13 Hz.) We 
introduced a parameter a in the exponential part because par- 
ticle acceleration models generally do not provide a robust 
prediction concerning the electron spectrum in the cutoff re- 
gion. We adopt a common value of a = 0.8; the choice has 
little influence on other parameters. Table[5]lists the spectral 
parameters that well describe the knot SEDs shown in Fig. [3] 
Also, in Fig.[6]we show the relative contribution of the second 
component, namely f 2u /{f\ v +hv), in the infrared to ultravi- 
olet part. It should be noted that the second component is 
dominant (> 80%) already in the R band flux and therefore 
the "optical" jet traces the high-energy radiation. For knot A, 
the first component has a luminosity of 1.3 x 10 43 erg s" 1 (in 
the frequency range 10 9 -10 14 Hz) without taking account of 
relativistic beaming. The minimum luminosity in the second 
component of knot Ais 3.4 x 10 43 ergs _I (10-10 18 Hz), and 
can be much higher if the spectrum extends well beyond the 
X-ray domain. 

The spectral index of the low-energy component is in the 
range of ct\ = 0.63-0.76, as determined from the radio spec- 
tra. Remarkably, the cutoff frequency is quite similar among 
the knots, v c = (4-6) x 10 13 Hz, which differs from previous 
modeling that de duced a progressi ve decrease of the cutoff 
frequencies (see Uester et al.l 12005. and references therein). 
The difference in the deduced cutoff frequency results from 
our inclusion of the IRAC fluxes into the modeling. The sec- 
ond spectral index for knots A and Bl is a 2 = 0.7-0.8; for 
knots CI, C2 and Dl, the second index is steeper, a 2 ~ 1 .0. It 
should be emphasized that the second index of knot A agrees 
well with the local spectral index determined by X-rays alone, 
in other words a 2 ~ a x . Interestingly, the steepness of the in- 
dex seems to increase with decreasing X-ray intensity of the 
knots. 



There are some caveats for the interpretation of the SEDs 
that we constructed. As mentioned earlier, knot B 1 in the op- 
tical shows a different morphology as compared to the radio 
one. Therefore, the first and second components should orig- 
inate in different emission volumes. Also, the X-ray flux as- 
signed for knot Bl could be contaminated by adjacent, bright 
knots. 

4.2. Physical Parameters for the Low-Energy Component 

The cutoff frequency of the low-energy synchrotron emis- 
sion that we determined in the previous section repre- 
sents the maximum energy of the associated electron pop- 
ulatio n. According to stand ard synchrotron theory (e.g., 
Ginzburg & Svrovatskii 1965), the maximum energy can be 
estimated as 



:0.28 



"1+z" 


0.5 


' " ) 




1 


v 0.1mGy 



-0.5 

( "* ) 

V10 14 Hz/ 



0.5 



TeV. 



(2) 

where 5 is the Doppler factor, z = 0.158 is the redshift, and 
B is the comoving magnetic field strength. The Doppler fac- 
tor is defined as S = [T (1 - f3 cos 6)]~ l with (3c the velocity 
of the jet, T = (1 - /3 2 )~ 1/2 the bulk Lorentz factor of the jet, 
and 9 the observing angle with respect to the jet direction. It 
is difficult to determine the Doppler factor, and therefore to 
what extent the observed jet radiation is enhanced by Doppler 
beaming. Furthermore, the jet may be decelerating along the 
observed jet portion resulting in a declining Doppler fa ctor 
( Georganopoulos & Kazanas 2004; Tavecc hio e t al. 2006). A 
likely range is 1 < 5 < 10 (see however 34. 4> . If equipar- 
tition between the radiating electrons and magnetic fields is 
assumed, without taking into account any relativistic motion 
(i.e., 5 = 1), one c an obtain an estimate of B = 0.1-0.2 mG 
(Jeste r et alJl2005l) . Therefore, given the cutoff frequencies 
of v c ~ 5 x 10 13 Hz, the maximum energy can be deduced as 
£max ~ 0.2 (fi-4<5)-° 5 TeV, wh ere B-a =B/(0.1 mG). Accord - 
ing to the prescription given in Harris & Krawczvnski ( 2002), 
the eq uipartition condition yields B cx <S -1 (see. lStawarz et alJ 
2003, for slightly different scaling). Then, an estimate of E m - dx 
is rather robust being nearly independent of both 6 and B. 

Synchrotron cooling in the source should affect the spec- 
tral shape provided that the cooling time is shorter than the 
formation timescale of the radiating knots. The synchrotron 
cooling time (in the jet comoving frame) for electrons with 
energy E = jm e c 2 emitting on average at a frequency v (ob- 
served) can be expressed as 



jm e c 

tsyn= (4/3)ca T t/ B 7 2 ' 

c n 0.5 



;2400 



(3) 



1+2 



B 



0.1 mG 



-1.5 



10 14 Hz 



-0.5 



yr,(4) 



where aj is the Thomson cross section, and Ub = B 2 /(8tv) is 
the energy density of the magnetic field. Equation (j2ji was 
used in relating E = ^m e c 2 with v, by replacing E m!ix with 
E and v c with (v/0.29). The source age, though quite un- 
certain, would be fo ~ 10 5 — 10 7 yr. Therefore, the observed 
synchrotron spectra from the radio to optical frequencies are 
likely to be formed in the "cooling regime", namely ? syn < fo, 
which means that a continuous injection of electrons with an 
acceleration spectrum of Q(E) oc E~ S[ (in the power-law part) 
over the timescale fo gives rise to the energy distribution of 
synchrotron-emitting electrons N(E) oc t syn Q(E) cx £'" (si+1) . In 
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FIG. 7. — Spitzer-HST-Chandra composite image of the jet of 3C 273. The colors are coded as follows: Spitzer "deconvolved" 3.6 pm (red), HST "UV excess" 
/o.3um -0.05 X fi.ft^m (green), Chandra 0.4-6 keV (blue). The VLA radio (2 cm) contours are superposed on the image, with the strongest radio source H2 being 
truncated. There are two distinct types of radiating knots; the inner knots show hard optical spectra (green) and strong X-rays (blue) while the outer knots show 
bright infrared (red). 



this case, according to the relation s\ = 2a\, the observed 
index of a.\ ~ 0.7 is translated into the electron accelera- 
tion index of s\ ~ 1.4, which is coincident with the hard- 
est possible index (s = 1.5) in the n onlinear shock accelera- 
tion theory (Malkov & Drurv 2001). We implicitly assumed 
partic le acceleration taking place in each knot (Jest er et all 
1200 ll) . where relativistic electrons are accumulated without 
escaping. A similar treatment was applied to the hotspot of 
3C 273 by Mei senheimer & Heavens! (11986) to incorporate a 
shock-acceleration theory into the emerging radiation spec- 
trum. There may be occasions for which the above argu- 
ment is not applicable; for example, in the case of significant 
beaming, the synchrotron cooling time would be much longer 
(f syn oc S 2 for the equipartition field) and the realization of the 
cooling regime would not be justified for radio-emitting elec- 
trons. 

4.3. Synchrotron Interpretation of the High-Energy 
Component 

For the inner knots (A and B) we have found that the opti- 
cal and X-ray fluxes can be modeled as a common power-law 
spectrum of spectral index — 0.7-0.8. Also, the optical to 
X-ray slope agrees with the X-ray spectral index, especially 
for knot A, ai ~ a x . One implication of this finding is that 
the optical and X-ray fluxes must be explained by the same 
emission mechanism. In 5j2] we mentioned that both the ra- 
dio and optical jet emission are linearly polarized to a s imilar 
degre e of ~ 15% llRoser & Meisenheime3 ll99ll iRoser et alJ 
1996), which has been taken as evidence for a synchrotron 
origin for the optical radiation in the inner knots. Now this in 
turn would imply a synchrotron origin for the X-ray radiation 
as well. Therefore we first discuss the optical-X-ray com- 
ponent in terms of synchrotron radiation produced by high- 
energy electrons. 

According to the prescription in 34.21 the power-law slope 
of the energy distribution of the accelerated electrons is S2 = 
2«2 in the synchrotron-cooling regime. Then the spectral in- 
dex of the high-energy component, 012 ~ 0.7-0.8 for the inner 
knots, yields S2 — 1.4—1.6. Interestingly, the index of acceler- 
ation is almost same between the low- and high-energy syn- 
chrotron components: S2 ~ Si (because of «2 ~ a\). On the 
other hand, the two indices differ in the outer knots: at2—l, 
corresponds to S2 — 2. 

The apparent difference in the cutoff frequency between 



the low- and high-energy components likely reflects a differ- 
ence in the maximum energy of electron acceleration. The 
cutoff frequency associated with the high-energy component 
is not constrained by the current data, but must be greater 
than 10 18 Hz (~ 4 keV). Consequently, from equation (j2ji, 
the lower limit to the maximum electron energy is E max > 
30 (B- 4 ST°- 5 TeV. 

The presence of two different cutoff energies suggests two 
acceleration modes that can be characterized by different ac- 
celeration rates in the knots. The acceleration timescale can 
be expressed as f acc = £r g /c, with a gyroradius r g = E/(eB). 
The rate of acceleration is parametrized by the factor £ > 1 
which is, as frequently assumed, taken to be energy indepen- 
dent. By equating f acc = f syn — namely, balancing the acceler- 
ation and synchrotron loss rates — one obtains the maximum 
attainable energy limited by synchrotron losses: 

£ma ^6xlO-r-(^)°' 5 eV. (5) 

Using equation this can be translated into the correspond- 
ing cutoff frequency, 

v c ~ 4.5 x 10 22 r 1 Hz. (6) 

1+z 

Therefore, for the second synchrotron component to extend 
beyond 10 18 Hz, the condition £ 2 < 4 x 10 4 should be fulfilled. 
On the other hand, for the low-energy synchrotron compo- 
nent, an extremely large value of £i ~ 10 9 can be inferred 
via a similar argument. Interestingly, a study of multiwave- 
length spectra of small-sca le blazar emission also gives large 
£ dlnoue & Takaharal 19961) . 

In the theory of diffusive shock acceleration, for rela- 
tivistic shocks, the factor £ corresponds roughly to the ra- 
tio of the energy density of regular magnetic fields to that 
of turbulent fields which resonantly scatter accelerating rel- 
ativistic particles. Then a different value of £ means a 
different development of turbulent magnetic fields near the 
shock. The inferred large value of £i , following a discus- 
sion in Biermann & Strittmatter (1987), may be understood 
in terms of steep turbulence spectra, e.g., of Kolmogorov- 
type I(k) oc k~ 5 / 3 [where I(k) is the magnetic energy density 
per wavenumber k in the turbulent field], though it should be 
noted that I(k) oc k~ l was silently assumed above to obtain 
energy-independent £ for simplicity. 
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There may be two distinct types of shocks characterized 
by different £ in the knots giving rise to their "double- 
synchrotron" nature. Alternatively, another acceleration 
mechanism may be operating, accounting in a natural way for 
a different acceleration rate for the second component. For 
instance, turbule nt acceleration m ay occur in the shear of jet 
boundary layers (Ostrowski 2000), generating the second syn- 
chrotron component (Stawarz & Ostrowski 2002). In highly 
turbulent shear layers, the acceleration rate can be estimated 
as £ ~ (3~^, where (3\ denotes the Alfven velocity in units of 
c. Provided that the jet is composed of ele ctron-pro ton plasma 
with the comoving number density n (ISikora & Madeiski 
2000), the Alfven velocity in the jet can be written as 



/3a = 



B 



= 7.3 x 10" 



B 



0.1 mG 7 Vl0- 4 cnr 



-0.5 



(7) 
(8) 



The condition £2 < 4 x 10 4 needed to produce the X-rays can 
be fulfilled if /3a > 5 x 10~ 3 , which seems reasonable for the 
typical jet parameters. 

An intriguing implication of interpreting the X-ray emis- 
sion as being due to a synchrotron process would be the pos- 
sibility that protons in the jet can be accelerated to very high 
energies. In the case of protons, synchrotron losses do not 
limit the acceleration process unless they reach ultra-high en- 
ergies. Instead, a condition set by f acc = f c , where t c is the light- 
crossing timescale of the knot along the jet (typically ~ 10 
kyr), would give a conservative estimate of the maximum at- 
tainable energy of protons. Here we use a light-crossing time 
rather than the "age" of the jet, probably to ~ 10 5 — 10 7 yr, to 
conservatively account for an escape loss which takes at least 
longer than a light-crossing time. This condition leads to 



;2.8 x 10 20 r' 



B 



0.1 mG / V 10 kyr 



(9) 

eV, (10) 



and consequently E p max > 0.7 x 10 16 eV for £ < 4 x 10 4 and 
typical values of B and t c . Remarkably, if £ ~ 10 (close to 
the so-called Bohm diffusion in the framework of diffusive 
shock acceleration), ultra-high-energy protons with energies 
ffn. max ~ 1Q 19 eV ca n be accelerated in the jet of 3C 273 
(see lAharonianl2002tlHonda & Hondal2004 for detailed the- 
oretical arguments). Even with a large value of £ (> 100), 
the shear acceleration would still be capable of generating 
£p.max ~ 10 19 eV, since the turbulent acceleration in the shear 
layer may be available for much longer timscales. The frac- 
tion of jet power that goes to high energy protons depends on 
the particle content of the jet and the details of the physics of 
particle acceleration, both of which are not well known. 

It is interesting to note that if protons are indeed acceler- 
ated in the jet to very high energies, £p. ma x ^ 10 18 eV, the 
second synchrotron component can be pr oduced directly by 
protons, i.e., proton synchrotron radiation ( Aharonian 2002), 
rather than by electrons. For the same particle energy, the 
characteristic frequency of proton synchrotron radiation, 2/ p>c , 
is less than that of electron synchrotron radiation by a factor 
of (m e /m p ) 3 . Therefore, by using equation 0, the spectrum 
of proton synchrotron has a cutoff frequency of 



'p,c 



:2.1 x 10 1 



1+z 



(- 

V mG 



10 18 eV 



which is located above the X-ray domain. Thus protons rather 
than a second population of electrons, may explain the sec- 
ond synchrotron component. Unlike the electron-synchrotron 
case, the diffuse nature of the emission (at least in the optical) 
can be understood more easily thanks to much longer cooling 
timescale for protons. A magnetic field strength of > 1 mG is 
needed in order to keep the power spent to accelerate protons 
at a reasonable level, and at the same time to render the accel- 
eration timescale shorter and the diffusion timescale longer 
jAharonianl2002l) . By adopting B = 10 mG, 6 = 1, an escape 
loss timescale of f esc = 1.4 x 10 7 (£ , p /10 14 eV)~ 1/ 2 yr, and a 
contin uous proton injection over to = 3 x 10 7 yr, lAharonianl 
(2002) obtained the acceleration power for protons of L p ~ 
10 46 erg s" 1 to explain the optical-X-ray spectrum of knot A 
in 3C 273 in terms of proton synchrotron radiation. 

4.4. Beamed IC Model for the High-Energy Component 

Previous explanations of the second, X-ray-dominated 
component of the SEP for jet knots have emphasized 
the be amed IC scenario (iTavecchio et alJl2000t ICelotti et alJ 



Hz, (11) 



2001). Specifically, relativistically beamed IC emission in the 
CMB photon fiel d has been invoked to m odel the X-ray emis- 
sion in 3C 273 (Sambrunaet al. 2001), as well as in other 
quasar jets. This model requires a highly relativistic bulk ve- 
locity of the jet, with a Doppler factor (5^10 out to nearly 
megaparsec distances. A direct connection between the opti- 
cal and X-ray fluxes in the inner knots A and Bl and also be- 
tween the ultraviolet and X-ray fluxes in the outer knots (see 
Fig. |3 led us to consider that the same mechanism should 
be responsible for the optical and X-ray emissions, forming a 
distinct high-energy component. If one adopts the beamed IC 
model for the X-rays, the optical fluxes in the inner knots and 
the ultraviolet fluxes in the outer knots should be explained 
also by the beamed IC radiation. This implies that the energy 
distribution of electrons has to continue down to very low en- 
ergies of E m i n ~ m e c 2 without a cutoff or break. 

In this scenario, one still should account for the optical po- 
larization. In the case of "cold" electrons in the jet, Comp- 
ton up-scattering off the CMB by an ultrarelativistic jet (so- 
called "bulk Comptonization") can yield in principle a high 
degree of polarizati on for the standard choice of 9 ~ T" 1 
( Beg elman & SikoralH987t> . However, for relativistic elec- 
trons with E ~ 2m e c 2 as appropriate for the optical IC emis- 
sion in this jet, the degree of polarization should be signifi- 
cantly reduced (Poutanen 1994). The present data of the ra- 
dio and optical polarization (Roser et al. 1996) may require a 
tuned set of physical conditions to reconcile with the beamed 
IC scenario. 

In this respect, a possible caveat would be that in addi- 
tion to the CMB, the synchrotron radiation of the jet itself 
might serve as seed photons for Compton up-scattering; if 
the highly relativistic jet accompanies slowly moving outer 
portions, the synchrotron radiation that comes from the slow- 
ing part can be relativistically enhanced in the fast moving 
jet, providing the effective photon fields for Comptoniza- 
tion ( Ghisellini, Tavecchio & Chiaberge 2005, as discussed 
for small-scale jets). If the observed radio flux is largely con- 
tributed by a slow, non-relativistic layer, the energy density of 
the synchrotron photons (as measured in the observer's frame) 
could reach ~ 10" 12 erg cm" 3 , being comparable to the CMB 
energy density. This possibility introduces a complication to 
the polarization argument. 

Within the framework of the beamed IC radiation, the 
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FIG. 8.— Broadband SEDs of the head region of the 3C 273 jet: H2 (filled 
circles) and HI (open circles). The photometric data are given in Table □ 
The curves represent the single-component model (see the text). 



physical parameters that can reprodu ce the SEDs of the 
3C 273 jet have been obtain ed (e.g.. iSambruna et ail 1200 It 
lHarris & Krawczvnskill2002l) . Compared to a typical value 
of E m j n ~ 20m e c 2 in the previous work (e.g., Sambruna et al. 
1200 ll) . a smaller minimum energy of the electron distribu- 
tions, say E m [ n ~ m e c 2 , is required by the present data. Note 
that this increases the estimate of the jet kinetic power by an 
order of magnitude if the jet power is dynamically dominated 
by cold protons. Let us model the SED of knot A with a s yn- 
chrotron plus beamed IC model of Tavecchio et al. ( 2000) as- 
suming 8 = r and a knot radius of 1 kpc (in the jet frame of 
reference). For a Doppler factor of 8 = 30 and magnetic field 
B = 5 /iG as roughly corresponding to an equipartition condi- 
tion, we obtain a jet kinetic power of Lj et ~ 1 x 10 48 erg s" 1 , 
if the number of the observed electrons is balanced by the 
number of cold protons in the jet. However, such a large 
Doppler factor of 8 = 30 seems problematic, since relatively 
large contributions from the un-beamed components are ob- 
served in the core emission, such as an ultraviolet bump (see 
ICourvoisierl 1 9981 and references therein), and the accretion- 
disk em ission c omparable t o the jet emission in the X-ray 
regime ( Grandi & Palumbo 2004). Instead, if we adopt a 
smaller Doppler factor of 8 = 10 (together with E m \ n = 3m e c 2 
to account for the optical emission), the energy density of 
electrons exceeds that of magnetic fields more than three or- 
ders of magnitude, and the jet kinetic power reaches an un- 
comfortable level of Ljet — 1 x 10 49 erg s" 1 . To summarize, 
the beamed IC model, in its simple form, does not easily find 
a satisfactory set of parameters, particularly for the case of 
knot A. 

Finally, the upper limit to the infrared fluxes from the quiet 
portion of the jet upstream of knot A is important to constrain 
the structure of the jet in the context of th e beamed IC model 
JGeorganopoulos et alJl2005t lUchivama et alJl2005l) . Unfor- 
tunately, we need additional IRAC observations to derive a 
reliable upper limit there because of the large systematic un- 
certainty in PSF subtraction in the region close to the quasar 
core. 

4.5. The Head Region 

In Fig. [8] we show the SEDs obtained for HI and H2 in 
the head region. Note that no significant X-ray emission was 
found in this region. We have modeled the SEDs using a 
power law with an exponential-type cutoff, namely by equa- 
tion with K2 = 0, as drawn in Fig. [8] The spectral parame- 
ters used for the model are indicated in Table|2] Interestingly, 
the spectrum of H2 shows a notable flattening at the optical 
and ultraviolet wavelengths, as first reported by iJester et afl 



(2005). 

The radio spectra are steeper than those of the upstream 
knots, with a ~ 0.9. It is likely that particle acceleration act- 
ing in the head region would be terminal shocks that have 
differen t properties as compared to upst ream shocks. As ar- 
gued by Meisenheimer & Heavens ( 1986), if the effect of syn- 
chrotron cooling is taken into account, the index seems con- 
sistent with an acceleration spectrum of Q oc E~ 2 , which is the 
standard spectrum of diff usive shock acceleration. Wit hin the 
picture proposed by Mei senheimer & Heavens! dl986l) . how- 
ever, the excess fluxes at optical and ultraviolet wavelengths 
seen in the spectrum of H2 cannot be easily understood. These 
excesses may represent a second synchrotron component just 
as in the case of the upstream knots. 

5. CONCLUSIONS 

We have presented results from our Spitzer IRAC obser- 
vation of the 3C 273 jet at wavelengths 3.6 and 5.8 /im, 
combined with the photometry with the VLA radio, HST op- 
tical, and Chandra X-ray data. Our multiwavelength anal- 
ysis led us to conclude that the flat optical emission in the 
X-ray-dominated knots originates in the high-energy power- 
law component, which also accounts for the X-ray emis- 
sion. The agreement between the optical-X-ray slope (02) 
and the X-ray spectral index (a x ) throughout the jet supports 
this picture. On the other hand, the radio to infrared spectra 
can be expressed by a power law with an exponential cut- 
off at ~ 5 x 10 13 Hz. The two distinct radiation components, 
namely low-energy (radio-infrared) and high-energy (optical- 
X-ray) emission, have similar power of ~ 10 44 erg s" 1 for the 
entire jet volume (without a beaming correction); the power in 
the second component becomes noticeably higher if its peak 
position is located far beyond the X-ray domain. The relative 
importance of the two components changes along the jet (see 
Fig-0- In the inner, X-ray-dominated knots, the high-energy 
component overwhelms the low-energy one. 

The second component can be attributed to either syn- 
chrotron radiation by a second population of high-energy 
electrons (or protons), or the beamed IC emission by the 
radio-emitting electrons. In the first case, the double syn- 
chrotron nature may arise from the presence of distinctively 
different acceleration processes (e.g., shock and turbulent 
acceleration). In the context of the origin of extragalactic 
cosmic-rays, it is interesting that a faster acceleration mech- 
anism producing the second component is capable of accel- 
erating cosmic-ray protons to energies 10 16 eV < E pmax < 
10 19 eV. On the Comptonization interpretation of the high- 
energy emission, we argued that the present polarization data 
seems problematic. Future sensitive measurements of polar- 
ization are quite important to draw firmer conclusions about 
the origin of the high-energy emission. 

The identification of the synchrotron component peaking at 
~5x 10 13 Hz in the X-ray-dominated knots of the 3C 273 jets 
calls for a follow-up deep observation with the Spitzer IRAC 
with better configurations. Our Spitzer cycle-3 proposal (PI: 
Uchiyama) requesting follow-up deep observations of 3C 273, 
as well as the observations of 10 other jets, has recently been 
approved. Ultimately, the James Webb Space Telescope, a fu- 
ture large infrared-optimized space telescope, will have an ex- 
cellent capability to explore the jet emission in the IRAC band 
with much improved resolution and sensitivity. 

An interesting implication of this work is that a num- 
ber of optical jets in radio-loud quasars discovered so far 
may be directly related to the X-ray emission rather than to 
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the radio emission, as in the 3C 273 jet. For example, in 
lUchivama et al.l fc005l) . we have associated the optical flux 
of the jet knot in the quasar PKS 0637-752 with the highest 
energy part of the synchrotron component. However, the op- 
tical flux could be attributed equally to the low energy part of 
the second component, given the result presented here. This 
issue should be addressed in future analyses based on mul- 
tiwavelength observations of C/ianofra-detected quasar jets. 
[In this context, we note that our discussions are not con- 
cerned with low-lumino sity FR I jets such as M 87 (e.g., 
iPerlman & Wilsonl 12005). in which the X-ray emission ap- 
pears to lie on the extrapolation of the emission at lower fre- 
quencies, and is widely believed to be of synchrotron origin.] 
It is clear from our results that the measurements over as wide 
a wavelength coverage as possible in the infrared-optical do- 
main are of great importance in studying the particle acceler- 
ation mechanisms in the large-scale jets of powerful quasars. 
Unfortunately, quasar jet emission in the optical region of the 
spectrum is very faint, and only a few optic al slopes for the 
powerful quasars ha ve been measured (Ridgway & Stockton 
1997; Cheung 2002). Further Spitzer observations are there- 



fore particularly helpful even though the spatial resolution is 
not optimal for the study of jet emissions. Also, a careful mea- 
surement of X-ray spectral index and its comparison with the 
optical fluxes is valuable to progress toward the understand- 
ing of the dominant X-ray emission mechanism operating in 
the quasar jets. 
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